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A model of neutrino oscillations is presented that has only three degrees of freedom and is consis- 
tent with existing data. The model is a subset of the renormalizable sector of the Standard-Model 
Extension (SME), and it offers an alternative to the standard three- neutrino massive model. All 
classes of neutrino data are described, including solar, reactor, atmospheric, and LSND oscillations. 
The disappearance of solar neutrinos is obtained without matter-enhanced oscillations. Quantita- 
tive predictions are offered for the ongoing MiniBooNE experiment and for the future experiments 
OscSNS, NOvA, and T2K. 
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I. INTRODUCTION 

Local Lorentz invariance is a basic feature of our best 
existing theory, which is the Standard Model (SM) cou- 
pled to General Relativity (GR). This theory is widely 
believed to be the effective low-energy limit of an under- 
lying structure that unifies quantum physics and grav- 
ity at the Planck scale, Mp ~ fO^^ GeV. Direct exper- 
imentation at this scale is infeasible, but sensitive mea- 
surements might detect suppressed low-energy signals of 
the expected new physics, including small violations of 
Lorentz and CPT symmetry [1, 2]. At presently attain- 
able energies, these unconventional effects can be char- 
acterized in the language of effective field theory [3] . Ex- 
tending the GR-coupled SM by adding all terms that in- 
volve operators for Lorentz violation and that are scalars 
under coordinate transformations results in a general ef- 
fective field theory called the Standard-Model Extension 
(SME). The leading terms in this theory include those 
of the SM and GR, together with ones violating Lorentz 
and CPT symmetry constructed from SM and GR fields 
[4, 5]. 

Neutrinos are natural quantum interferometers, and so 
neutrino oscillations represent sensitive phenomena with 
which to conduct searches for physics beyond the min- 
imal SM, including Lorentz and CPT violation. Com- 
pelling experimental evidence for oscillations exists, and 
they are conventionally attributed to neutrino masses. In 
the standard three- neutrino massive model [6] , the effec- 
tive hamiltonian for neutrino propagation and oscillation 
is a 3x3 matrix determined by six parameters, consist- 
ing of two squared- mass differences AttIq, Am^^j^^, three 
mixing angles 9i2, 023, ^13, and a phase S. This model 
can reproduce the observed features of the solar-neutrino 
suppression [7-11], the atmospheric- neutrino oscillations 
[12-15], and the KamLAND results [16]. However, the 
model cannot reproduce the LSND signal [17]. 

The presence of Lorentz and CPT violation in the 
underlying theory can produce additional or alternative 
sources of neutrino oscillations in the low-energy effec- 
tive Lagrange density, which are contained in the neu- 
trino sector of the SME. Even if attention is restricted 
to three generations of active neutrinos and dominant 
effects, the SME involves many terms and many types 



of effects [18]. These oscillation effects are controlled 
by coefficients for Lorentz violation forming various 3x3 
matrices denoted (a^)"^, (cl)"^, etc., where a, (3 are 
Lorentz indices and the matrix indices a, b label neu- 
trino flavor and span e, fi, r. It is elegant to regard these 
coefficients as arising from spontaneous Lorentz violation 
[19], which is consistent with Riemann geometry [5] and 
may be ubiquitous in effective field theories [20]. Some 
of the neutrino coefficients for Lorentz violation act to 
mimic aspects of conventional mass-induced oscillations, 
while others predict qualitatively new features of oscil- 
lations such as sidereal variations, direction dependence, 
unconventional spectral behavior, Lorentz-violating see- 
saws, neutrino- antineutrino mixing, and more. 

The LSND collaboration recently analyzed the 
sidereal-time dependence of their reported 9^ appearance 
signal [21], using the short-baseline approximation [22]. 
Since this signal cannot be incorporated within the stan- 
dard three-neutrino massive model, a possible explana- 
tion is Lorentz violation. The scale required to gener- 
ate a neutrino-oscillation signal at the LSND energies 
E ~ 10-2 GeV was found to be of order IQ-^^ GeV 
for (ai)g^ and [cl)"^E. This result is compatible with 
the dimensionless ratio ~ 10"^^ of electroweak to Planck 
scales that might be expected to control suppression of 
Planck-scale physics. The sensitivity achieved is compa- 
rable to other searches for Lorentz and CPT violation us- 
ing interferometric techniques in high-energy physics in- 
volving K oscillations [23] , D oscillations [24] , and Bd os- 
cillations [25-28] , including ones with sidereal variations 
[29]. More generally, the SME sensitivities attainable 
with various types of neutrino searches can be compara- 
ble [18] to those achieved in the numerous experiments 
in astrophysics, atomic physics, optical physics, nuclear 
physics, and particle physics [1, 2]. 

Each coefficient for Lorentz violation leads to an en- 
ergy behavior different from that of a neutrino mass, so 
estimates for the size of possible Lorentz violation may 
be found by considering the energy dependence of the 
existing neutrino-oscillation data [30-34]. With this ap- 
proach, profitable experiments for certain types of coeffi- 
cients would involve ultra-high-energy neutrinos [35] be- 
cause the long-baselines and high-energies involved would 
yield high sensitivity to certain small Lorentz-violating 
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terms. However, it has been shown that some spe- 
cial combinations of coefficients for Lorentz violation can 
mimic a mass-like energy dependence through a Lorentz- 
violating seesaw mechanism, even if the hamiltonian con- 
tains no mass terms. The so-called 'bicycle' model [36] 
offers a simple example with only two coefficients for 
Lorentz violation. Despite having no neutrino mass in 
the hamiltonian, the model reproduces many features of 
the neutrino-oscillation data through a pseudomass term 
generated by a seesaw. 

In this work, we describe a model for neutrino oscil- 
lations based on the SME that, like the bicycle model, 
uses one CPT-odd coefficient and one CPT-even coeffi- 
cient. However, we also introduce a mass term, so the 
hamiltonian has three degrees of freedom. This 'tandem' 
model has two Lorentz-violating seesaw mechanisms that 
interplay in interesting ways in certain energy regimes. 
We demonstrate that the model reproduces the features 
of all reported neutrino-oscillation data, including the 
LSND signal. The model and some basic theoretical top- 
ics are discussed in Sec. H, while the application of the 
model in the context of neutrino-oscillation data is pre- 
sented in Sec. IH. Wc find that observed L/E behaviors 
are reproduced for all baselines L and energies E, and 
the solar-neutrino result is obtained without the use of 
matter-enhanced oscillations [37]. 



II. TANDEM MODEL 

In this section, we first discuss some criteria that are 
useful in guiding the construction of neutrino models in 
the presence of Lorentz and CPT violation. The tandem 
model is then obtained, and some of its basic properties 
are presented. 



A. Criteria 

In the standard three- neutrino massive model [6], the 
effective 3x3 hamiltonian describing the propagation and 
oscillation of neutrinos of energy E takes the form 

{hSMjab = Edab H ^ • (Ij 

The squared-mass matrix (m|]yj)a6 can be specified by 
two squared-mass differences, three mixing angles, and a 
phase. Except for the LSND signal, the current neutrino- 
oscillation data are consistent with two large mixing an- 
gles and one approximately zero. It follows that (tosm)^^ 
can be written phenomenologically in the four-parameter 
matrix form 

/O \ 

^^Im ~ Uiu Am| UsM (2) 



with 

/ 1 ^ \f cos^ia sin^ia 0\ 

C^SM = cos 6*23 sin 6*23 -sin 6*12 cos 6*12 
\0 -sin6i23 COS6I23 / \ 1/ 

(3) 

The existing data are consistent with the parameter val- 
ues Ami - 8-0 X 10-5 eV2, Am^^^ ~ 2.5 x IQ-^ eV^, 
6I12 34°, and 6*23 45°. 

In this paper, we develop an alternative model that 
adequately describes all the existing neutrino-oscillation 
data, including the LSND signal. The effective hamil- 
tonian for neutrino propagation in this model contains 
a term involving conventional neutrino mass, together 
with an admixture of Lorentz-violating terms. To date, 
no compelling evidence for Lorentz violation exists, so 
a model of this type is of immediate interest only if 
it is more attractive than the conventional picture in 
some other respects. The standard three-neutrino mas- 
sive model (1) has a solid foundation in renormalizable 
quantum field theory, and it is consistent with all data 
other than LSND using only four parameters. Moreover, 
its mass scales < 0.1 cV arc compatible with a seesaw 
origin [38, 39]. Therefore, if Lorentz violation is to be 
invoked, it is desirable to consider models that (i) are 
based on quantum field theory, (ii) involve only renor- 
malizable terms, (iii) offer an acceptable description of 
the neutrino-oscillation data, (iv) have any mass scales 
< 0.1 eV for seesaw compatibility, (v) involve fewer pa- 
rameters than the four used in the standard picture, (vi) 
have coefficients for Lorentz violation consistent with a 
Planck-scale suppression < 10"^'', and (vii) can accom- 
modate the LSND signal. 

How challenging is it to satisfy these seven criteria? 
To satisfy (i) it suffices to focus on the SME, since this 
provides a general field-theoretic framework for Lorentz 
violation. Satisfying (ii) requires restricting attention 
to Lorentz-violating operators of dimension four or less. 
However, this restriction is highly nontrivial when taken 
in conjunction with (iii). Coefficients for Lorentz viola- 
tion satisfying (ii) generate neutrino oscillations that are 
either independent of E or proportional to E, in sharp 
contrast to the observed behavior proportional to 1/E 
required to satisfy (iii). 

One potential solution to this problem is illustrated in 
the bicycle model [36], which is based on the minimal 
SME and describes well the behavior of solar and atmo- 
spheric neutrinos [40] through a Lorentz-violating see- 
saw mechanism. The bicycle model has no mass terms 
and only two coefficients for Lorentz violation suppressed 
by the Planck scale, so (iv) is irrelevant and (v), (vi) 
are satisfied. The observed E~^ dependence of oscilla- 
tions at large energies emerges as a combination of the 
Lorentz-violating E^ and E^ dependences. However, at 
lower energies the model predicts a direction-dependent 
constant-energy signal, which may be excluded by Kam- 
LAND data. Also, the LSND signal remains unexplained. 
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B. Hamiltonian 

The goal of the present work is to provide an exphcit 

example of a model satisfjdng all seven criteria (i)-(vii). 
To satisfy (i) and (ii), we adopt Lorentz- violating terms 
from the minimal SME but omit for simplicity the pro- 
vision for neutrino- antincutrino mixing. In this context, 
the effective hamiltonian for neutrino propagation takes 
the form [18] 



With this assumption, the neutrino effective hamilto- 
nian can be written in the form 



{Ks)ab = E5ab + 



2E 



+^[{aLrp^-{cLrp^PuU. (4) 

Since the coefficients [abjab associated with CPT-odd 
operators in the Lagrange density, the effective hamil- 
tonian (h'^Q)ab for antineutrino propagation is obtained 
by reversing the sign of the coefficients {aLj^b- The co- 
efficients (ol)^5 have dimensions of mass, while {chjab 
are dimensionlcss. In obtaining Eq. (4), possible gravi- 
tational couplings [5, 41] have been disregarded, and the 
coefficients (az,)^^ and [chlab assumed to be space- 
time constants. If these coefficients originate in sponta- 
neous Lorentz violation, they have companion Nambu- 
Goldstone fluctuations that could be interpreted as the 
photon [42], the graviton [43], or additional neutrino in- 
teractions [44], but these effects are secondary in the 
present context and are disregarded in this work. 

The occurrence of the neutrino three- momentum p in 
the hamiltonian (4) means that the oscillation physics 
in the chosen inertial frame typically depends on the di- 
rection of neutrino propagation [18]. The number of de- 
grees of freedom can be significantly reduced if this com- 
plication is avoided. One possibility is to suppose the 
model is rotationally invariant. In a Lorentz- violating 
theory, this requirement can be implemented only in a 
single special inertial frame. If this frame is identified 
with that of the cosmic microwave background radia- 
tion, for example, then direction-dependent effects are 
still present in neutrino-oscillation experiments because 
the solar system is moving with respect to this frame. 
However, the rotational, orbital, and translational mo- 
tions of the Earth are essentially nonrelativistic in the 
approximately inertial Sun-centered frame [45] relevant 
for neutrino experiments, so in practice any direction- 
dependent effects in models of this type are suppressed 
to parts in a thousand or more. Another possibility is to 
suppose the model does have large direction-dependent 
effects but that they are irrelevant for describing the 
available experimental data. The point is that reported 
neutrino-oscillation results hiq typically obtained by inte- 
grating data taken over long time periods, so the average 
values of the direction-dependent coefficients may suffice 
even in theories with large direction-dependent effects. 
In what follows, we disregard direction-dependent effects 
as a reasonable first approximation, without precluding 
their existence in a more detailed treatment. 



{Ks)ab 



Edab + + K)a6 - l{CL)abE. (5) 



In practice, terms proportional to the unit matrix can be 
disregarded because they produce no oscillation effects. 
Note, however, that mass terms proportional to the unit 
matrix may play a role in ensuring stability and causality 
of the underlying theory [46]. 

To maintain the possibility of satisfying criteria (iii) 
and (vii) while keeping the number of degrees of freedom 
small, we further restrict attention to effective hamiltoni- 
ans of the form (5) adapted from the simple scheme of the 
bicycle model [36]. In that model, the high-energy pseu- 
domass is created from a Lorentz- violating seesaw mech- 
anism involving a (ci)-type coefficient in the on-diagonal 
(^cfi)ee component together with off-diagonal (ai)-type 
coefficients. To incorporate also the observed KamLAND 
L/E dependence [16] and generate effects that can repro- 
duce the LSND signal [17], a second 'tandem' Lorentz- 
violating seesaw can be introduced that operates at low 
energies. This can be triggered by adding another on- 
diagonal entry that is located in the {hes)TT component 
and involves an (m^)-type parameter. These considera- 
tions suggest limiting attention to the special case of the 
hamiltonian (5) taking the form 



' 3 (cL)ee 
{aL)ne 



Jer 



(6) 



Since this hamiltonian is hermitian and the elements are 
real, it is symmetric. The model (6) therefore has five 
degrees of freedom. 

To satisfy criterion (v), we must reduce the degrees of 
freedom to fewer than four. We find that the current 
data may be described with the simplifying assumptions 
that all the (aL)-type coefficients arc identical, (ai)e^ = 
{<1l)ht = {0'L)Te- For simplicity, following the suggestive 
notation of Ref. [18], we write mrr = m, ql = a, and 
— 4(cL)ee/3 = c. The effective hamiltonian for neutrinos 
becomes 



"TM 




(7) 



This is the tandem model, which depends on only three 
independent degrees of freedom and is the focus of the 
remainder of this work. Note that the presence of CPT 
violation implies that the corresponding effective hamil- 
tonian for antineutrinos is 



cE —a —a 
—a —a 
-a -a 'm?/2E 



(8) 



In Sec. Ill, we show that a choice for the three de- 
grees of freedom can be made that respects criteria (iv) 
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Feature 


standard model [oj 


Bicycle [36] 


Tandem 


Lagrange-density formulation compatible with other physics 


yes 


yes 


yes 


renormalizable terms only 


yes 


yes 


yes 


generations of active neutrino species 


3 


3 


3 


number of sterile neutrinos 











degrees of freedom in effective liamiltonian 


4 to 6" 


2 


3 


global description of neutrino-oscillation data except LSND 


yes 


yes" 


yes 


description oi LSND result 


no 


no 


yes 


matter-enhanced oscillations needed 


yes 


yes 


no 


seesaw-naturaJ masses (it any) < 0.1 eV 


yes 


none 


yes 


coefficients for Lorentz violation (if any) < 10~^^ suppressed 


none 


yes 


yes 


CPT-odd operators 


no 


yes 


yes 


large direction-dependent effects 


no 


yes 


no= 


neutrino-antineutrino mixing'* 


no 


no 


no 



"Depending on whether Bis and S are included. 
''Except perhaps KamLAND (cf. Sec. II A). 
"^In the isotropic approximation. 
'*See Ref. [18] for models with this feature. 



TABLE I: Some attributes of the standard three-neutrino massive model, the bicycle model, and the tandem model. 



and (v) and that yields a global description of the cur- 
rent neutrino-oscillation data inclnding the LSND signal. 
This means that the tandem model satisfies all seven cri- 
teria (i)-(vii), making it a useful alternative candidate 
model for neutrino oscillations. In fact, the structure of 
neutrino oscillations in the tandem model is remarkably 
rich despite the presence of only three degrees of freedom. 
This is a consequence of the double Lorentz- violating see- 
saw and the accompanying CPT violation, and it pre- 
dicts a variety of observable phenomena in future exper- 
iments. Table I provides a summary comparison of some 
attributes of the standard three-neutrino massive model, 
the bicycle model, and the tandem model. 



C. Properties 

To analyse neutrino mixing in the tandem model, we 
diagonalize the hamiltonian (7) using a 3x3 unitary mix- 
ing matrix U, 

h^u = U'iEtmU, (9) 

where Etm is a 3 x 3 diagonal matrix containing the h^y^ 
eigenvalues Ej, J = 1,2,3. The description of neutrino 
oscillations then depends on the mixing matrix elements 
Uja and the eigenvalue differences Ajk = Ej — Ek- 

Since /iJ^m is nondiagonal and includes terms with dis- 
tinct dependences on the neutrino energy E, both Uja 
and Ajk have intricate energy behavior. Moreover, the 
E dependences of the corresponding quantities for the 



antineutrino hamiltonian (8) are different. As a result, 

the oscillation probabilities for neutrinos and antincutri- 
nos can vary strongly and distinctly with E. For ex- 
ample, we show in the next section that the dependence 
of Uja on E results in energy variations for solar pp- 
and ^B-neutrino oscillations analogous to those arising 
via matter-enhanced effects [37] in the standard three- 
neutrino massive model. Similarly, the dependence of 
Ajk with E produces behavior matching the data for 
both the KamLAND reactor and the Super-Kamiokande 
atmospheric neutrinos. Also, an oscillation signal for 
LSND emerges. 

To perform the diagonalization explicitly, we used two 
different procedures: numerical matrix diagonalization, 
and analytical solution of the cubic eigenvalue equation. 
The two procedures yield the same results. For the nu- 
merical work, we adopted the simple Jacobi method [47] . 
An appropriate orthogonal matrix (Jacobi rotation) is 
created and applied to the hamiltonian to eliminate one 
off-diagonal component. This is repeated in turn for the 
other two off-diagonal elements, thereby completing one 
Jacobi sweep. Five sweeps yield diagonalization at suffi- 
cient precision, permitting calculation of the eigenvalue 
differences Ajk and mixing-matrix elements Uja- For 
the analytical method, the standard del Ferro-Cardano 
solution of the cubic equation was used to obtain the 
three eigenvalues and the corresponding eigenvectors and 
to construct the mixing matrix. Some details of this so- 
lution are given in the Appendix. In practice, the diag- 
onalization procedure was repeated for each value of the 
neutrino energy E of interest. 
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The neutrino-oscillation probability in the SME frame- 
work is derived in Ref. [18]. For the tandem model, all 
coefBcients in /i^j^ are real and the oscillation probability 
for neutrinos reduces to the simple form 

P^^^^, = - 4 ^ UjaUjbUKaUxb Sin^ {^AjkL), 
J>K 

(10) 

where L is the baseline distance. The oscillation proba- 
bility Pjpj^jj^ for antineutrinos is given by an expression 
of the same form, but with Uja and I^jk obtained by 
diagonalizing h'i^j^ instead. 

Since the effective hamiltonian /i^j^ is symmetric, 
Pjy^^^i^ = P^^^^^, as may also be seen from Eq. (10). 
This is a property associated with T invariance. In fact, 
the Lorentz violation arising via nonzero isotropic coef- 
ficients a and c preserves T symmetry. The presence 
of the coefficient a implies CPT violation, but detecting 
this violation via P^^^^,^ ^ Pp^^p^ may be challenging in 
certain energy regimes, as discussed in the next section. 
Since the tandem model is CPT violating but T invari- 
ant, it breaks CP symmetry. More generally, Lorentz 
violation maintaining rotation invariance and controlled 
by coefficients of the and Cl types preserves T sym- 
metry but violates CP through a combination of C and 
P breaking [18, 48]. Note that this origin of CP violation 
is qualitatively different from that in the standard three- 
neutrino massive model, which arises via a phase in the 
mixing matrix U [6]. 

Finally, we offer a few remarks about model build- 
ing. At present, no completely satisfactory theory of 
Icpton and quark masses is available. In the standard 
three-neutrino massive model, understanding the neu- 
trino mass matrix offers some unique challenges associ- 
ated with the origin, values, and stability of the masses 
and mixing angles. The standard picture provides par- 
tial answers by invoking one or more seesaw mechanisms 
in a grand-unified theory, perhaps with supersymmetry 
[49, 50]. In this context, provided criterion (iv) is satisi- 
fied in the selection of the value of rh?, the tandem model 
is on a roughly comparable footing to the standard three- 
neutrino massive model. However, the required structure 
of the neutrino mass matrix is somewhat different, which 
offers interesting scope for model building. The essential 
issue for the tandem model is to explain the dominance 
of the component (m^)^^ in the mass matrix in Eq. (4). 
This can arise naturally in some cases. For example, in a 
simple SO(IO) grand-unified theory the neutrino Dirac- 
mass matrix can be proportional to the quark mass ma- 
trix [50, 51], so invoking a seesaw mechanism can produce 
neutrino masses proportional to the square of the quark 
masses. Under these circumstances, the large mass of 
the t quark can ensure that the neutrino mass matrix 
is dominated by (m^)^^, as is effectively assumed in the 
tandem model. Analogous model-building issues exist for 
the coefficients for Lorentz violation c and a in the tan- 
dem model, and these represent an interesting open area 
for future work. 



III. APPLICATION 

We find that the tandem model provides a good de- 
scription of existing neutrino-oscillation data with the 
following values for the mass parameter and the two co- 
efficients for Lorentz violation in Eqs. (7) and (8): 

im^ = 5.2 X 10-3 eV^ 
a = -2.4 X lO"^'' GeV, 
c = 3.4 X 10"^^. (11) 

Note that this choice of m respects the cosmological con- 
straint on neutrino masses [52] . Also, the values for a and 
c are consistent with the results extracted from LSND 
data [21, 22]. 

The choice of values for m^, o, and c is guided by 
the presence of the double Lorentz-violating seesaw in 
the tandem model. Taken independently, each seesaw 
generates distinct asymptotic effects, but by continuity 
the two must merge in some energy regime. Near the 
merger scale, the neutrino behavior is controlled by the 
interplay of both seesaws, and interesting effects appear. 
The values (11) are chosen to yield these merger effects 
aromid 10 MeV, while still satisfying criteria (iv) and 
(vi). It is possible that other values for m, a, and c exist 
that are compatible with the neutrino-oscillation data. 

In this section, we begin by describing some general 
features predicted by the model with the values (11). We 
then offer a comparison to existing results from a variety 
of neutrino-oscillation experiments. 



A. Features 

Some general features of the tandem model with the 
choices (11) are shown as a function of energy for neutri- 
nos in Fig. 1 and for antineutrinos in Fig. 2. These fig- 
ures reveal a relatively complicated energy dependence, 
arising because the energy dependences of the eigenvalue 
differences Ajk and the elements of the mixing matrix 
Uja are more involved than the simple power law of the 
standard three-neutrino massive model. As can be seen 
by comparing Figs. 1 and 2, the oscillation lengths and 
probabilities of neutrinos and antineutrinos differ in de- 
tail. This occurs because the coefficient a is associated 
with CPT violation and enters with the opposite sign in 
the two cases. 

Figures 1(a) and 2(a) display curves in the logL-logi? 
plane that establish experimental sensitivities to oscilla- 
tions involving particular flavors. The oscillation length 
L for each type of oscillation is determined by the asso- 
ciated eigenvalue difference Ajk, and the corresponding 
curve is the solution to the equation AjkL = tt that fixes 
the occurrence of the first oscillation maximum. An ex- 
periment is sensitive to oscillations of a particular type 
if it lies in a region of the logi-logi^ plane above the 
corresponding curve. 
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FIG. 1: Neutrino behavior in the tandem model, (a) Loca- 
tion of the first oscillation maximum in the log L-log E plane. 
The dashed, dotted, and dash-dotted lines correspond to A21, 
A32, and A31, respectively. For comparison, the thin and 
thick solid lines show the location for Aniatm and AmQ, re- 
spectively, in the standard three-neutrino massive model, (b) 
Values of —4UjaUjbUKaUKb versus energy for « oscil- 
lations. The dashed, dotted, and dash-dotted lines correspond 
to {J,K) = (2, 1), (3,2), and (3, 1), respectively, (c) Same for 
Vij. <-> Vt oscillations, (d) Same for z/e <-> Vt oscillations. 



In the standard thrcc-ncutrino massive model, the os- 
cillations depend on L/ E, and the condition for the first 
oscillation maximum is described by a straight line with 
logL oc logE in the logL-logi? plane. In Lorentz- 
violating models, a coefficient of the type also pro- 
duces a straight line for the the first oscillation max- 
imum, but with zero slope: logL is constant. A coeffi- 
cient of the Cl type produces a straight line with negative 
slope, logi oc — logi?. Neither of the latter two behav- 
iors is observed in the neutrino-oscillation data. However, 
Lorentz-violating models with both types of coefficients 
can produce complicated energy dependences via seesaw 
effects [18]. These effects produce the sensitivity curves 
for the tandem model shown in Figs. 1(a) and 2(a). 

Note that the tandem model does exhibit a limiting 
straight-line sensitivity behavior with logL oc logE for 
A21 at atmospheric-neutrino energies and baseline dis- 
tances. This featiirc is also present in the bicycle model, 
and it matches the atmospheric-neutrino data [36, 40]. 
At high energies, the limiting behavior of both A32 and 
A31 = A21 + A32 is also linear, but with logi oc — logii^. 
In the low-energy limit, the sensitivity behaves as logL oc 
— logE for A21 and as logi oc logE for A32 and A31. 

Figures 1(b)- (d) and 2(b)-(d) show the energy depen- 
dence of the quantities —WjaUjbUKaUKb in Eq- (10) 
that control the amplitude of the oscillation probabili- 
ties. Unlike the standard three-neutrino massive model, 
these quantities contribute to the energy dependence. It 
turns out that they permit a description of the energy 
dependence in the solar-neutrino data without resorting 
to matter-enhanced oscillations [37]. 




log E(MeV) log E(MeV) 



FIG. 2: Antinoutrino behavior in the tandem model. The 
descriptions of (a)-(d) match those of Fig. 1. In (a), the dotted 
and dash-dotted lines appear to coincide but in fact remain 
distinct for all as explained in the Appendix. 



B. Comparison to data 

Some generic predictions of the tandem model for 
typical experimental sensitivities are displayed in Fig. 
3. These plots present information about appearance 
and disappearance probabilities in the \ogL-\ogE plane 
for both neutrinos and antineutrinos. They provide an 
overview of the general predictions of the model and re- 
late them to some existing and future experiments. Some 
details for the various classes of experiments are given in 
the following subsections. 

In Figs. 3(a) and 3(b), regions above the curves corre- 
spond to disappearance probabilities greater than 10%. 
In Figs. 3(c) and 3(d), regions above the curves cor- 
respond to appearance probabilities greater than 0.1%. 
The figures also contain rectangles representing the re- 
gions of sensitivity of particular experiments. The tan- 
dem model predicts a signal in an experiment for oscilla- 
tions of a particular type when its rectangle overlaps the 
region above the oscillation curve. 
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vacuum solar 
MINOfRg 





n 

Bugey, □ v^^v^ 
Gosgen q ^ 




r □ v'^v^i 

r d) □ v^v - 



1 2 3 4- 
log E(MeV) 



1 



12 3 4 
log E(MeV) 




ECMeV) 



FIG. 3: Oscillation probabilities in the logi-logiJ plane for 
(a) V disappearance, (b) v disappearance, (c) v appearance, 
and (d) v appearance. In (a) and (b), the displayed regions 
have disappearance probability P > 10%. Regions for dis- 
appearance of different flavors are distinguished by different 
hatchings. In (c) and (d), the displayed regions have ap- 
pearance probability P > 0.1%, and the regions for different 
flavors are again distinguished by different hatchings. Ap- 
proximate L-E sensitivity ranges for various oscillation ex- 
periments are displayed as labeled rectangles. In (a), the 
rectangle for solar neutrinos lies well above the plot, near 
logL ~ 11. Short-baseline experiments with S > 1 GeV are 
denoted SB. 



1. Solar 

In the standard three-neutrino massive model, vac- 
uum solar-neutrino oscillations show no energy depen- 
dence because AttIq ^ 10"^" eV^ and the mixing ma- 
trix is energy independent. In the tandem model, how- 
evcir. the long-baseline limit of vacuum oscillations ex- 
hibits an energy dependence arising from the mixing ma- 
trix Uja, despite the loss of energy dependence in the 
factor sin^(AjifL/2). 

The Ve survival probability on the energy scale rele- 
vant to the solar neutrino problem is shown in Fig. 4. 
This model yields an oscillation probability of 50% for 
pp neutrinos, which have a continuous energy spectrum 
with endpoint at 0.420 MeV. The oscdllation probability 
is also 50% for ''Be neutrinos, which generate two mo- 
noenergetic lines with 90% at 0.862 MeV and 10% at 
0.384 MeV. For the neutrinos, which have a contin- 
uous energy spectrum with endpoint at 15.04 MeV, the 
model yields a probability of ~ 40%. These values are 
in reasonable agreement with the energy dependence ob- 
served in the existing solar-oscillation data. 



FIG. 4: Solar-neutrino oscillation probabilities in the tandem 
model, (a) Curves representing survival probabilities for Ue 
(solid), Ue — » (dashed), z/^ — » Ur (dotted), and Vt — > I'e 
(dash-dotted), (b) Survival probability of Ue in the tandem 
model (solid line) and in the standard three-neutrino massive 
model with a basic matter-induced effect (dashed). Note the 
different energy scales. The effects of experimental position 
and energy resolution are not shown. 



In the tandem model, matter-enhanced oscillations [37] 
play essentially no role in the description of solar neutri- 
nos. In general, matter-enhanced solar-neutrino oscilla- 
tions can be viewed as arising from an effective coeffi- 
cient for Lorentz violation of the form [18] (aL,eff)ee = 
V^Gpne, where He is the number density of electrons in 
the Sun. An analytical approximation for rig [53] yields 
(aL,eff)ee - 1-87 X io-i7e-io-54fl/flo MeV, where R is 
the distance from the core and Rq is the solar radius. 
In the tandem model, substantial effects can arise only if 
cE < (aL,eff)ee, Or £; < 0.56" ^O-^"*^/^© MeV. For much 
of the solar-neutrino spectrum, this energy is too small 
to allow appreciable effects on the oscillation probability. 
This result is confirmed by numerical calculations. 

We remark in passing that an interesting feature of 
the tandem model in the 10 MeV region is a reduced 
suppression of the Ug probability relative to that of Vg- 
This originates in CPT violation via the a coefficient, and 
it may be relevant for r-process nucleosynthesis following 
core-collapse in supernovae [54] . 



2. Atmospheric 

The atmospheric-neutrino results are well described by 
the tandem model. At high neutrino energies, the effec- 
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L/E (m/MeV) 



FIG. 5: Survival probabilities for atmospheric neutrinos as a 
function oi L/E for (a) i'^ and (b) Curves arc shown for 
the tandem model (solid) and for the standard three-neutrino 
massive model (dashed) . The effects of experimental position 
and energy resolution are not shown. 
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L/E (m/MeV) 



FIG. 6: Survival probabilities for long-basclinc reactor an- 
tineutrinos. (a) Survival probabilities as a function of E for 
i>e (solid), i/e (dashed), —> Vr (dotted), and Vt — » 

(dash-dotted), (b) Survival probabilities for as a function 
oiLjEm. the tandem model (solid) and in the standard three- 
neutrino massive model (dashed) . The effects of experimental 
position and energy resolution are not shown. 



tive tandem hamiltonian for neutrinos becomes approxi- 
mately 



TM 



/ cE 
\ 



a 

a a J 



a 

a a 



(12) 



with a similar simplification for /Ixm- The c and a co- 
efficients combine to produce a Lorentz- violating seesaw 
mechanism [18], in which the model effectively reduces to 
a two- flavor limit with a pseudomass term [36] . The only 
significant probabilities are for the transitions i/^ Ut 
and v^j, Vr, which involve A21 and are large only at en- 
ergies above about 100 MeV. All amplitudes other than 
U2ij.U2tUiij,Uit approach zero at high energies. 

A comparison of atmospheric-oscillation probabilities 
in the tandem model with those in the standard three- 
neutrino massive model is provided in Fig. 5. In both 
models, the disappearance probability is maximum in the 
region with L/E = 300-600 m/MeV, which corresponds 
to the results reported by the Super-Kamiokande exper- 
iment [14]. This feature is common to both u^-v^ and 
j^/i-i'T oscillations. The tandem model contains CPT vio- 
lation, but it is invisible in the high-energy region. These 
results are also consistent with the recent atmospheric 
measurement from the MINOS experiment [15]. 



3. Long-baseline reactor 

At low neutrino energies, the effective tandem hamil- 
tonian for antineutrinos becomes approximately 



/ 



■m?/2E I 



(13) 



The component rr? jlE creates another Lorentz- violating 

seesaw [18]. Here, it results in large ^ f'^ oscillations 
involving A21, with no contribution from other antineu- 
trino flavors. The only oscillations are proportional to 

The oscillatory shape reported by the KamLAND ex- 
periment [16] is reproduced by the choice of values (11). 
Figure 6 shows the S-g survival probability as a function 
of energy and LjE. 

Unlike the situation for atmospheric neutrinos, the 
CPT violation in the tandem model emerges in this L-E 
range as a difference between and survival proba- 
bilities. However, this difference is unobservable in the 
KamLAND experiment because the sources produce only 
v^. Note that the oscillation probability is sub- 

stantial only at low energies. This means the tandem 
model predicts a null signal for future long-baseline v^- 
appearance experiments such as NOvA [55] and T2K [56] . 
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Experiment 


Ref. 


Oscillation channel 


Source 


Baseline 


Energy 


Status 


Bugey 


[57] 

L J 




,-> Y 




reactor 


15,40 m 


3 MeV 


null result 


Gosgen 


[58] 




ff, — > A 




reactor 


38, 46, 65 m 


3 MeV 


null result 


Palo Verde 


[59] 








reactor 


750, 890 m 


3 MeV 


null result 


CHOOZ 


[60] 




. Y 

— > y\- 




reactor 


1 km 


3 MoV 


null result 


KARMEN 


[61] 




Vp, i'e 




accelerator 


18 m 


40 MeV 


null result 


LSND 


[17] 

L J 




l^H — ^ 




accelerator 


30 m 


40 MeV 


signal 


OscSNS 


[62] 




l^H — ^ 




accelerator 


60 m 


40 MeV 


planning 


MiniBooNE 


[63] 

L J 




I'ti — » l^e 




accelerator 


550 m 


800 MeV 


ongoing 


CDHS 


[64] 

L J 




Z/ji — » v\ 




accelerator 


130 m 


1 GeV 


null result 


BNL-E776 


[65] 


Vu 


^ fe, i-'p — » 


Ve 


accelerator 


1 km 


1.4 GeV 


null result 


CHORUS 


[66] 




ff^ Vr 




accelerator 


600 m 


27 GgV 


null result 


NOMAD 


[67] 




Vr 




accelerator 


600 m 


45 GeV 


null result 


CCFR 


[68] 




Ve, 


Ve 


accelerator 


1 km 


140 GeV 


null result 


NuTeV 


[69] 




-> Ve, V^ 


Ve 


accelerator 


1 km 


150 GeV 


null result 



TABLE II: Recent and near-future experiments involving short-baseline oscillations at reactors and accelerators. The notation 
Va — > Vb indicates an appearance search, while Va ^ X indicates a disappearance search. 

4- Short-baseline reactor and accelerator 

In recent years, many experiments with short baselines 
[17, 57-69] have searched for oscillations using neutrinos 
from reactor and accelerator sources. The results arc 
summarized in Table II. The predictions of the tandem 
model relevant at these baselines and energies are dis- 
cussed in this subsection. 

In the tandem model, oscillation probabilities at short 
baselines are typically small, so they are undetectable 
in reactor-disappearance experiments with sensitivities 
of order 10%. This can be seen from Fig. 3. Conceivably, 
a sensitive reactor experiment with a baseline of > 1 km 
might observe a signal predicted by the tandem model in 
the highest-energy bins of the v^ spectrum. 

For short baselines < 1 km, oscillations of neutrinos at 
high energies > 100 MeV are strongly suppressed, as can 
also be seen from Fig. 3. The oscillation behavior in the 
high-energy region is similar to that of atmospheric neu- 
trinos, and any experiments with high-energy neutrinos 
that are insensitive to atmospheric-parameter oscillations 
are also insensitive to oscillations in the tandem model. 

In contrast, for short-baseline low-energy accelerator 
experiments, the nonzero mass term in the tandem model 
results in a small oscillation probability that lies on the 
edge of experimental sensitivity. The i/^ v^ and 
v^ Ve oscillation probabilities for KARMEN, LSND, 
OscSNS, and MiniBooNE arc shown in Fig. 7. The tan- 
dem model yields an oscillation probability on the order 
of 0.05-1% in the detectable energy range for all these 
experiments. 

In particular, the model yields an acceptable value for 



T -1 - a) 
ilO ^ 

^ -2' 

10 r 



-3- 



10 



^.10 

T -1 ^ c) 

^10 r 

^ -2 = 
10 r 

-3- 

10 ^ 



10 




b) 





1 10 10^ 10^ 
E (MeV) 



1 10 10^ 10^ 
E (MeV) 



FIG. 7: Oscillation probabilities as a function of E for neu- 
trinos (solid) and antincutrinos (dashed) in (a) KARMEN, 
(b) LSND, (c) the proposed OscSNS experiment, and (d) the 
currently running MiniBooNE experiment. The effects of ex- 
perimental position and energy resolution are not shown. 
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the oscillation probability of (0.264± 0.067± 0.045)% re- 
ported by LSND. The oscillation probability is smaller 
for the KARMEN experiment due to the shorter base- 
line, and it lies below the experimental sensitivity. The 
tandem model predicts substantial appearance signals for 
Ue and De in OscSNS and MiniBooNE. Note that these 
signals fall off rapidly with increasing energy. 



the entries are all real and given by ft-ee = cE, = 0, 
Htt = rr? I2E, and /le/* = ^re = ^^it = 

The first step is to solve the cubic eigenvalue equation. 
This takes the form 

+ oA^ + 6A + c = 0, (A2) 

where 



IV. SUMMARY 

In this work, we have constructed a model of neu- 
trino oscillations involving only three degrees of freedom 
that describes all features of existing oscillation experi- 
ments, including the LSND result. This tandem model 
uses renormalizable terms in the SME, with one mass pa- 
rameter and two coefficients for Lorentz violation. The 
value of the mass parameter is compatible with a conven- 
tional seesaw mechanism, and the coefficients for Lorentz 
violation have Planck-suppressed magnitudes that are 
consistent with the measurements recently reported by 
LSND. The observed solar-neutrino suppression and en- 
ergy dependence are achieved without employing matter- 
enhanced oscillations. 

The tandem model makes quantitative predictions that 
can be tested in the near future with results from the 
MiniBooNE experiment and from various planned short- 
and long-baseline appearance experiments such as Osc- 
SNS, NOvA, and T2K. There may also be subsidiary sig- 
nals involving sidereal variations or other unconventional 
physics. In any case, the tandem model represents a can- 
didate alternative to the standard three-neutrino massive 
model, and its existence reconfirms the key role of neu- 
trino experiments in the ongoing exploration of physics 
beyond the Standard Model. 
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APPENDIX A: DIAGONALIZATION 

The general symmetric effective hamiltonian repre- 
sented by the 3x3 matrix 




(Al) 



can be analytically diagonalized, /leff = If^EeeU, by solv- 
ing the cubic eigenvalue equation for the eigenvalues Ej 
and then obtaining the diagonalizing matrix Uja- In the 
effective hamiltonians (7) and (8) for the tandem model. 



a = -tr(/icfr), 

b = i[tr(/ieff)]2-itr(/i2g), 

c = -det(/ieff). (A3) 

For the tandem model, these quantities are 

a = -cE-m^/2E, 
b = cmV2-3a^ 

c = (cE T 2a -h rh'^/2E) . (A4) 

Note that b is energy independent and that CPT violation 
contributes only to the determinant quantity c. 

It is convenient to introduce the combinations of a, b, 
and c given by 

Q = i(a2-36), 

R = ^{2a^ -9ab + 27c), 

f R \ 



e 



cos 



(A5) 



In terms of these combinations, the energy eigenvalues 
Ej can be expressed as 



El 



cosl- 



E2 = -2^/Qcos(^ 



3«' 



+ 2Tr 



^3 = -2VQcos' ^ 1- 



(A6) 



In general, the solution to the cubic (A2) involves three 
unequal real eigenvalues Ej provided the discriminant 
D = R? — Q^ ]s negative definite. For the tandem model, 
the values (11) chosen for the model satisfy this condi- 
tion for all neutrino energies E. With these values, the 
discriminant D increases with i? to a maximum near 10 
MeV, then decreases as E increases. The appearance 
of the maximum around 10 MeV reflects the merger of 
the two Lorentz-violating seesaw mechanisms near that 
energy. Also, since the eigenvalues Ej of the tandem- 
model effective hamiltonian are distinct for all energies 
E, no two eigenvalue differences Aj^ coincide exactly 
for any given E. This fact is reflected in Figs. 1(a) and 
2(a): if two eigenvalue differences were to coincide ex- 
actly then the logarithm of the third would diverge, a 
behavior absent from the flgures. 
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The mixing matrix Uja is found to take the form 



and the normalization for each fixed J is 



U = 



where 



/ BiCi/Ni B2C2/N2 B3C3/N3 ^ 

CiAJNi C2A2/N2 CsAs/Ns 
\ AiBi/Ni A2B2/N2 AaBa/Ns J 



Aj = hij,r{hee — Ej) — hreheix, 
Bj = hre{hf_i^ — Ej) — hefihfj^r, 
Cj = henihrr — Ej) — h,j,rhre, 



(A7) 



(A8) 



Nj 



A^jBj + B'jC'j + CjAy 



(A9) 
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